###########################################################################
# Correlations of Authoritarian Items with Support, p.26

c.b <- cor.test(data$behave, data$liberal, use="complete.obs")
c.s <- cor.test(data$safety, data$liberal, use="complete.obs")
c.t <- cor.test(data$told, data$liberal, use="complete.obs")

auth <- (data$behave + data$safety + data$told)/3
c.m <- cor.test(auth, data$liberal, use="complete.obs")

cor.auth <- rbind(c(c.b$estimate, c.b$p.value),
       c(c.s$estimate, c.s$p.value),
       c(c.t$estimate, c.t$p.value),
       c(c.m$estimate, c.m$p.value))

cor.auth <- round(cor.auth, 2)
colnames(cor.auth) <- c("Corr", "p-value") 
rownames(cor.auth) <- c("behave", "safety", "told", "average") 

c.alpha <- round(cronbach.alpha(cbind(data$behave, data$safety, data$told), na.rm=T)$alpha, 2)
###########################################################################
